d x 2_ y 2 superconductivity in a generalized Hubbard model. 
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We consider an extended Hubbard model with nearest-neighbor correlated hopping and next 
nearest-neighbor hopping t' obtained as an effective model for cuprate superconductors. Using 
a generalized Hartree-Fock BCS approximation, we find that for high enough t' and doping, antifer- 
romagnetism is destroyed and the system exhibits d-wave superconductivity. Near optimal doping 
we consider the effect of antiferromagnetic spin fluctuations on the normal self-energy using a phe- 
nomenological susceptibility. The resulting superconducting critical temperature as a function of 
doping is in good agreement with experiment. 



I. INTRODUCTION 

One of the most challenging properties that any effec- 
tive one-band model for the high superconducting critical 
temperature (high-T c ) superconducting materials should 
explain is the experimentally observed dependence of T c 
with doping (x). Extended t — J models seem to be suc- 
cessful to interpret several important properties of these 
materials |] §. They were derived as effective Hamilto- 
nians in the strong-coupling limit (very large U) |I] , 
while in generalized Hubbard models obtained by simi- 
lar derivations, U is of the order of the unperturbed band 
width 8t ~ 3eV While numerical studies in t— J- like 
models unambiguously indicate that they exhibit sizable 
superconducting correlations with d-wave symmetry, in 
agreement with experiments H, this is not the case of 
the usual Hubbard model [jioj. Apart from quantitative 
reasons concerning the strength of the local Coulomb re- 
pulsion, it is of interest to study an accurate enough effec- 
tive model amenable to weak coupling many-body treat- 
ments. In particular, the observed pseudogap behavior 
for temperatures T c < T < T* in angle-resolved pho- 
toemission spectroscopy (ARPES) experiments jllj has 
been interpreted as a precursor effect of the antiferro- 
magnetic (AF) jl2| as well as the superconducting state 
p3| , H between other possible scenarios [|l5|,^6). For the 
attractive (negative U) Hubbard model, a good deal of 
research helped to elucidate the nature of the supercon- 
ducting transition and most of this work is based on 
one appealing feature of this model: its explicit attrac- 
tive interaction leading to a superconducting state in the 
Bardeen-Cooper-Schrieffer (BCS) approximation. Thus, 
even though the symmetry of the BCS gap is s-wave, 
this model has been used for a qualitative understanding 
of ARPES experiments |13|. Due to the lack of a gen- 
uine microscopic d-wave superconducting analogous of 
the of the negative-?/ Hubbard model, some phcnomcno- 



logical models with BCS-like interactions were studied 
@H|l9||. While the model proposed in Ref. Jlf] is built 
up by pure phenomenological interactions with several 
free parameters, in order to fit the experimentally ob- 
served T c vs x curve, the AF and van Hove picture [l9| ] 
has been proposed on the basis of a simplification of the 
t — J model in which holes are constrained to move only 
in one sublattice (up or down) in a Neel spin background. 
Nevertheless, this is an oversimplification for realistic 
J < t p0[ . The assumption of pairing coexisting with 
long-range AF order is shared with the strong-coupling 
polaron picture ||, which also reproduces qualitatively 
the observed T c vs x curve. However, nuclear magnetic 
resonance experiments indicate that the coherence length 
of the AF correlations £ is of only a few lattice sites in 
the optimally doped materials pl| . An effective attrac- 
tive separable potential has also been considered to study 
the d-wave superconducting pseudogap evolution p^ j . In 
this case, however, the AF correlations are not taken into 
account at all, while they play a relevant role for under- 
doped and optimally doped materials fl2| , |22|| . 

In this work, we calculate T c vs x of an effective 
one-band model for the high T c cuprates using weak- 
coupling techniques. We obtain beyond a certain doping, 
a stable paramagnetic phase with <i-wave superconduc- 
tivity (DWS). The mechanism for superconductivity is 
the electron-hole symmetric correlated hopping studied 
previously in one and two dimensions (2D) How- 
ever, previously in 2D, the next-nearest-neighbor hop- 
ping t' was neglected and only s-wave superconductivity 
was found for doping high enough to inhibit long-range 
antiferromagnetism. In simplified terms, this can be un- 
derstood as follows: on-site Coulomb repulsion U inhibits 
on-site pairs, but nearest-neighbor singlet pairs are fa- 
vored by the correlated hopping. In mean field, the de- 
pendence of this effective attraction in reciprocal space is 
proportional to (cosfca; — cosk y ) 2 for d symmetry and to 
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(cos k x +cos ky) 2 for extended s symmetry (see Section III 
or Ref . [pi| ) . This implies that <i-wave superconductivity 
is favored with respect to the s-wave one, for doping x at 
which the points X ((0,7r) and equivalent) lye near the 
Fermi surface. However, when t' — 0, this happens for 
x ~ and for these dopings, due to the perfect nesting 
of the Fermi surface, the antiferromagnetic instability is 
the dominant one and inhibits superconductivity. When 
t' is included, the perfect nesting is destroyed and the X 
points lye at the Fermi surface for finite doping, stabiliz- 
ing the d-wave solution. 

In Section II we briefly describe the derivation of the 
effective model and explain the meaning of the different 
terms. In Section III we calculate T c vs x using a general- 
ized Hartree-Fock (HF) BCS decoupling. Since the above 
mentioned decoupling exaggerates the range of stability 
of the AF phase, in Section IV we include the effect of 
the AF spin fluctuations (AFSF) for dopings at which 
long-range AF order is destroyed, by evaluating the nor- 
mal self-energy in the one-loop approximation with the 
phenomenological susceptibility of Ref. f2]}| . Section V 
contains a discussion. 

II. THE EFFECTIVE MODEL 

The model used is obtained from the mapping of the 
low-energy part of the Hilbert space of the three-band 
Hubbard model H 3 b As first suggested by Zhang 

and Rice [||, and later confirmed by analytical and nu- 
merical work of several groups ^-fl], when insulating 
cuprates are doped, the added holes form local singlets 
involving a Cu 3d x 2_ y 2 hole at a given site i (with cre- 
ation operator d\ a ) and a hole in the linear combina- 
tion a\ a of the four nearest-neighbor oxygen 2p CT orbitals 
with x 2 — y 2 symmetry. This Zhang-Rice singlet has the 

form [4(a| T 4; - a%4{) + Bd h d li + Ca l-\ a li\\°) and is 
mapped into the vacuum |0) at site i in the effective one- 
band model H. Similarly, the ground state of for 
one hole at cell i, (which has the form (Dd ia + Ea\ a )\Q)) 
is mapped onto ct_ CT |0) of H, where cJ CT is an effective 
electron creation operator. The vacuum of at cell i 
(which corresponds to full 3d and 2p shells) is mapped 
onto c^c|j0) in H. Calculating the matrix elements in 
the reduced Hilbert space, and retaining only the most 
important terms, the following effective Hamiltonian re- 
sults: 

H = U^n l] n ii - H^2{n^ + nn) 

i i 

~ X] ( C i<x C Jff + h.c){t AA (1 - n ia )(l - Tlja) 

<ij>a 

+ t B B n ia n jrT + tAB [«-«t(1 - n ja ) 

+ n ja (l - n i(7 )]} - t' c L c jV, (1) 

<ij'>cr 



where < ij > (< ij >) denotes nearest-neighbor (next- 
nearest-neighbor) positions of the lattice. Note that U 
is not directly related with a Coulomb repulsion, but 
represents the cost in energy of constructing a Zhang- 
Rice singlet from two singly occupied cells. It is lower 
than the difference between p and d on-site energies of 
H 3 b (usually called charge-transfer energy A). tAA rep- 
resents the hopping of a Zhang-Rice singlet to a singly 
occupied nearest-neighbor cell. The terms with ampli- 
tude tAB correspond to destruction of a Zhang-Rice sin- 
glet and a nearest-neighbor cell without holes, creating 
two singly occupied cells and vice versa. tsB describes 
the movement of an isolated hole in H^. Clearly, the 
amplitude of these three correlated hopping processes 
should differ in general. The dependence of the next- 
nearest-neighbor hopping on the occupation of the sites 
involved in the hopping is neglected. Note that in the 
simplest strong-coupling derivation leading to the t — J 
model Q, tAA amounts to the only kinetic term in the 
model, the one with tsB is mapped out of the relevant 
Hilbert-space while the one with tAB is treated in second 
order of perturbation theory to define the exchange term 
. More systematic strong-coupling derivations contain 
other terms as well as higher-order corrections ■ 

For t' = 0, the occurrence of metal-insulator transi- 
tions and superconductivity have been investigated for 
particular cases of this model [^3|,^4) , but a stable DWS 
phase has not been found before. We concentrate on the 
parameter regime tAB > tAA = tBB = t, which seems 
more adequate for the cuprates [pp|j2^|. The effect of 
t' Q is crucial to stabilize the rf-wave phase, because it 
removes the perfect nesting at half-filling and shifts the 
energy of the van Hove singularity (VHS) in the unper- 
turbed density of states (originated by the saddle points 
at X in the dispersion relation) away from the Fermi en- 
ergy at half filling. 

III. MEAN-FIELD APPROACH. 

The correlated hopping terms of the Hamiltonian (|l|) 
can be separated in one- two- and three-body contribu- 
tions, with coefficients t = tAA, £2 = tAA — tAB and 
<3 = 2tAB — tAA ~tBB respectively |m| - We treat Eq. 
(|]) within the generalized HF BCS approximation [p4[ . 
The term in t% contributes to the BCS solution in the 
d-wave as well as in the s-wave channels. The self- 
consistent parameters considered in the decoupling are 
(n iT ) - {riiy) = me lQ R *, r = (c| +(5(T c l(T ), ip = (c^c^) and 

<fis = (c\ + si c li), with S = x,y being vectors connecting 
nearest neighbors and (p x = ±tp y (^4). n = n^+n^ = 1—x 
is the particle density, L is the number of lattice-sites, 
Q = (ir, 7r) and indicates the lattice-position. The 
problem is reduced to a one-particle one with three pos- 
sibilities for the symmetry breaking perturbation: (i) AF 
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spin density wave (SDW) (with m/0 and tp — fx = 0), 
(ii) DWS (with m — tp — and <p x — — tp y ^ 0), and (iii) 
extended s-wave superconductivity (SWS) (with m = 0, 
tp ^ and </? x = if y 7^ 0). The possibility of coexistence 
of SDW and superconductivity is left out here, since a 
previous study indicated that a sizeable superconducting 
gap is not possible within our model in presence of long- 
range antiferromagnetism (^6|. For the three cases, the 
rcnormalized dispersion relation, effective hopping and 
effective chemical potential can be written in the form: 

£fe = — 2t e ff(cos k x + cos k y ) — At 1 cos k x cos k y , 
teff = t- t 2 n + t 3 [3r 2 + tp 2 + ipl - [n 2 - to 2 )/4] 
Hefj = fi- [Un/2 + U 2 T + At 3 [rn + 



*P(fx+fy)])- 



(2) 



The SDW, d-wave, and s-wave BCS order parameters 
are: 



A BCS 
^Dk 



A 



BCS 
Sk 



("2 +4i 3 T) m 

4t^Tip x (cos k x — cosfc y ), 
(8i 3 r - U)<ip - A{2t 2 + nts)tp a 
+ [4t 3 rtp x - 2(2t 2 + nt 3 )')p] 
(cos k x + cos k y ), 



first instability (as T is lowered) of the Hessian matrix 
formed by the second derivatives of f2 with respect to the 
independent variables [p8[. 



0.06 



0.04 



0.02 



0.00 




(3) 



respectively. For a given wave vector k, they coincide 
with half of the corresponding energy gap. For both, d- 
and s-wave superconducting solutions, the dependence 
of T c vs x is obtained from the linearized BCS gap equa- 
tions. For the d-wave case, T c is given by: 



f 



- — -TT tann — 
(2tt) 2 K 2T' 



4i 3 T (cos k x — cos k y ) 
2lfc 



cos fc T 



(4) 



where Ek = — f^effj &ud T is the temperature. For the 
s-wave solution, T c is the temperature at which \ ma x 
equals one, being \ m ax the largest eigenvalue of the ma- 
trix 




tanh(|^ 
tanh(^) 



C tfe COS fc-E 

2e k 



Jl0tanh(#)^ 
J^ tanh( ^)&-i 



where 



8i 3 r — U — 2{2t 2 + nt 3 )(cos k x + cos k y ), 
-4(2t 2 + nt 3 ) + 4t 3 r(cos fc 2 + cos fc„), (6) 



are the coefficients of tp and ^ in Ag£ s . This method 
of obtaining T c in second-order phase transitions when 
the thermodynamic potential Q depends on more than 
one parameter (here tp and f x ), has been used before 
p7f, and is equivalent to the usual one of finding the 



FIG. 1. d-wave BCS- critical temperature T c in units of t 
as a function of the doping x (for x < 0.5), for tAB = 2t (open 
symbols) and t ab = 1.5t (filled symbols) . Circles, squares and 
triangles correspond to t' = 0, — 0.2t, — 0.45£, respectively. 
The solid symbols for x > 0.5 show T c of the s-wave solution 
for tAB = 1.5 and the same values of i' as those considered for 
the d-wave. The dot (dot-dashed) line indicates the boundary 
of the SDW at T = for U = it (U = 6t), £ab = 1.5t 
(tAB = 2t) and £' = -0.45t 

In Fig. 1, we show the superconducting critical tem- 
perature as a function of doping x = I — n, where n is 
the number of electrons per site for both, d-wave and 
s-wave solutions. Most of the pairing terms, as well as 
the additional contribution to the usual U/2 term in the 
expression of the SDW gap, can be generated from the 
mean-field decoupling of an effective two-body interac- 
tion of the form 



(7) 



<i]> 



Paring in both d-wave and extended s-wave channels, is 
originated by the spin- flip terms of V e f t . For the s-wave 
case, there are some additional contributions, as can be 
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observed in (||) . From mere inspection of the three possi- 
ble order parameters (^) , it can be seen that A SDW and 
Ag^ s depend on U in such a way that the antiferromag- 
netic (s-wave-BCS) solution is enhanced (weakened) as 
U increases, whereas does not depend on U. Thus, 

the doping region where the d-wave BCS solution exists 
is limited only by the difference in energy with the two 
other competing instabilities of the Fermi liquid. As men- 
tioned in Section I, due to the wave-vector dependence 
of the superconducting order parameters, the shape of 
the Fermi-surface plays an important role in defining the 
doping regions at which each of the two BCS solutions 
are possible. In Ref. J24|] we have shown that for t' = 
and assuming a constant density of states, SWS exists 
for high enough x and low enough U, while DWS exists 
near half-filling. Within this region, however, the SDW 
solution has lower energy than the DWS one, being the 
difference between the energy of both solutions small for 
U ~ (as can be inferred from the form of the gaps, 
Eqs (||)). When t' is turned on, the VHS is displaced, 
lying at the Fermi energy for a finite doping. The incre- 
ment in the density of states at finite doping enhances 
the critical temperature as well as the doping range asso- 
ciated to the s-wave solution. However, as a consequence 
of the vector dependence of the effective interaction, the 
optimal doping for SWS is always high. In Fig. 1, it is 
shown that for tAB = 1-5 and U = it, the doping region 
with SWS is x > 0.6, i.e. beyond the range of doping ac- 
cessible experimentally in the cuprates. For the <i-wave 
solution, the maximum of T c occurs when the chemical 
potential coincides with the energy of the VHS. Due to 
the loss of perfect nesting and the closing of the indirect 
SDW gap for this doping, the SDW is weakened. Thus, 
even for large U , for sizeable t' there exists a doping 
range for which the d-wave BCS solution is stable and 
robust. We indicate in the figure the doping region for 
the existence of the SDW solution (m ^ 0) at T = 
for the particular values U — At (U = 6t), tAB = l-5i 
(tAB = 2t) and t' — — 0.45t. The Neel temperature at 
which m becomes zero is not calculated here since: it is 
determined by spin fluctuations, and the mean-field value 
is too large. For nearly realistic bare parameters, the 
doping range for an stable DWS (beyond the dot-dashed 
line in Fig. 1) corresponds to the overdoped regime of 
the cuprates, in agreement with experiments on Bi2212 
Pt which indicate that the superconducting state takes 
place in an otherwise ordinary Fermi liquid for overdoped 
compounds. 

In the optimally doped regime, spin fluctuations not 
only affect directly the SDW order parameter, destroying 
the AF long-range order, but also affect indirectly the su- 
perconducting order parameters. Their effect should be 
included in a more realistic calculation of the supercon- 
ducting critical temperature. This is the goal of the next 
section. 



IV. THE EFFECT OF SPIN FLUCTUATIONS 

Long-range antiferromagnetism is not expected to sur- 
vive within the whole range of doping predicted by the 
mean-field approach. For doping larger than x ~ 0.05 — 
0.1, the spin fluctuations are important and a Fermi liq- 
uid with strong AFSF is a more appropriate picture for 
this regime |^l| , p5[ . AFSF of any Hubbard-like micro- 
scopic model as the one considered here can be charac- 
terized by a spin susceptibility of the form 

x(Q,o) 



x(q,^) = 



(8) 



i + (q - Q) 2 e - iu/usf ' 

where the parameters scale with the AF correlation 
length £ as x(Q,0) ~ £ 2 , uj sf = r £ 2 , with r ~ 40meV 
pl| , p9[ . Assuming a phenomenological spin-fermion cou- 
pling, AFSF has been suggested to mediate the pairing 
in the cuprates ^5|. This picture is supported by the 
observed correlation between £ and T c in several super- 
conducting cuprates. Within the weak-coupling formal- 
ism, DWS due to correlated hopping is inhibited in a 
background with long-range AF correlations plfl . How- 
ever, the effective interaction Eq. (Q) provides an explicit 
channel for the coupling with collective AFSF within the 
doping region without AF long-range order, which has 
the same form as the phenomenological coupling used in 
Ref. HI. Thus, it might be expected that the AFSF 
would renormalizc the bare value of tAB ( < 1.5t for re- 
alistic values of the three-band parameters |J ) to higher 
ones [[22|]2"5|| . In what follows, we investigate how the 
BCS- T c vs x dependence is modified by the effect of the 
AFSF for the d-wave solution. 

For finite t', the Fermi surface contains hot spots (for 
which €k = Cfe+g). Fermions located in the neighborhood 
of these points are the most affected by AF correlations 
and exhibit a peculiar T-dcpcndcncc in the onc-particlc 
spectral properties, which is mainly determined by the 
magnitude of £/£ t h> and w 3 f/T |22p^], with £ th = v F /T, 
being vp the Fermi velocity. Hot spots are located at 
near X ((0, 7r) and symmetry related points), i.e. near 
the antinodes of the DWS gap and with energies close to 
the VHS (for which vp ~ 0). As T c vs x predicted by the 
BCS approximation, as well as the value of the maximum 
T c itself, depend on the behavior of the density of states, 
AFSF are expected to play some further important role 
apart from the eventual renormalization of the effective 
pairing interaction. The self-energy obtained from (@) is 



S(k, iu n ) = T 



E 



d 2 q 



s(q)<?(q) 



x(q,w m ) = - 



(2nr 

x(q,w m )G°(k + q, iuj n 
u ° duj Imx(q, u>) 



7T iv„ 



(9) 



where [G (k, iu> n )] 1 = iui n — e^, g(q) is an effective 
interaction between fermions and spin fluctuations, ujq 
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is a frequency cut off, v m — 2rmrT and ui n = (2n + 
l)nT. As usual, an effective coupling constant g' ~ 
5( c l)5( ( l)x(Qj 0)/C 2 is defined, which in the present case 
should be proportional to (U + At^r ) ^(|. The ensu- 
ing spectral function is -A(k, uj) = — ImG(k, u>)/tt, with 
[G(k,w)] _1 = u> — ik — S(k, uj). To examine how the 
changes in the behavior of A(k, u>) affect the T c vs x de- 
pendence, we consider the effect of AFSF using the BCS 
form of the anomalous self energy Af^ s given by Eq. 
(|J) and calculating the normal self energy in the one- 
loop approximation ( Eq. (|9|)). The resulting linearized 
gap equation, 

^ = tJ ^K C k S cosfe^e^+GCk,^) 

G(-k,-iw n ), (10) 
can be cast in real frequency as 

<e* = -j 4^2 A ?S S cos*! 

J -f- ZJ. 

where n f (u) = 1/(1 + cxp((w - [i ef )/T)). Eq. @ 
reduces to the linearized version of the usual BCS Eq. 
(Q) when A(k, w) = S(ui — ejt). As discussed in previ- 
ous works, assuming £ approximately constant or with 
a weak T- dependence, two different regimes due to the 
AFSF can be distinguished as a function of T in the be- 
havior of A(k,w) HIJUl: (i) For T « u> s f, quantum 
contributions dominate the behavior of the self-energy 
(^J) and j4(k, uj) exhibits Fermi liquid-like quasiparticle 
peaks, even for Fermi points near k^ s . Within this 
regime, no relevant qualitative changes in comparison 
with the BCS description are expected in the solution 
of the gap Eq. 0). (ii) For T » uj sf , and for k- 
points satisfying £ >> £th, classical effects (introduced 
by the m = O-Matsubara frequency in (^)) dominate. 
The AFSF can be considered as quasistatic and A(k} ls , uj) 
exhibits a shadow-band structure for large enough values 
of g' |^,^9|. This implies a transfer of spectral weight 
from low to high frequencies, and we expect an effective 
blurring of the large density of states near the VHS with 
a concomitant decrease of T c . 

The above qualitative issues are confirmed by our nu- 
merical calculations, as illustrated in Figure 2. We re- 
strict ourselves to the case with finite t' = — 0.45i, which 
reproduces the observed Fermi surface of YBCO and 
B2212 and to the case of tAB = 2t, for which the op- 
timal doping and the maximum T c predicted by the BCS 
approximation are in good agreement with experiments 
(for these parameters, the ratio t'/t e ff ~ 0.27). To com- 
pute T c we evaluate /x e / / and r from 

n=l + T J ^^ e ^ 0+ ReG(k, lw „), 




0.00 1 ' ' ' 1 

0.00 0.10 0.20 0.30 0.40 0.50 
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FIG. 2. Critical temperatures T c in units of t as a 
function of the doping x, for tAB = 2t and t' = — 0.45t, 
considering AFSF. The dashed line (crosses) corresponds 
to the BCS-solution in the thermodynamic limit (using 
fc-discretization) . Circles correspond to g' = 3t, and £ = 2.5 
(open ones) and £ = 10 (filled ones). Squares correspond to 
g' = Wt and £ = 2.5 (open ones) and £ = 4 (filled ones). 

and afterwards check if there exists (or not) a solution 
of the linearized gap Eq. (|i"(i|). A cut-off tt m = 432w s / 
was chosen for the i/ m -Matsubara summation in @ while 
another one Q n = 46i, assuming t — 250meV, was used 
for the w„-Matsubara summation in ( |To| ) and jl^). The 
sum over the tails was approximated by an Euler-Mc Lau- 
rin formula. The fc-integral in Eq. (|^) is calculated with 
a relative precision of 10 -6 , while a fixed finite mesh was 
used to evaluate the ^-integrals in Eqs. ([n]) and jl^). 
In the latter step, some precision is lost and the BCS- 
T c vs x, calculated in the thermodynamic limit, is not 
exactly recovered as indicated in Fig. 2. We considered 
£ = 2.5, (oj sf = 0.0256i) and £ = 4, (uj sf = O.Olf), 
which are supposed to be representative of YBCO and 
LaSrCuO, respectively pl| near optimal doping, and 
£ = 10, (u) s t = 0.0016t) which could be realistic only 
for strongly underdoped materials. Taking g' — 3t a 
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weak effect on T c is observed. For g' = lOt, which is 
expected to be representative of the cuprates [ p5|]29[ , a 
decrease in T c is observed as £ increases. For the case of 
£ = 10, g' = lOt (not shown in the figure) the maximum 
T c is below 0.012*. 

V. DISCUSSION 
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We have shown that for high enough next-nearest- 
neighbor hopping t', d-wave superconductivity is stabi- 
lized for doping x ~ 0.35 at the mean-field level in the 
correlated hopping model (|l|). The corresponding crit- 
ical temperature T c has the right order of magnitude 
(~ 70 — lOOif). Our calculations are in 2D and this 
result may be, of course, affected by the fluctuations of 
the superconducting state, which has not been taken into 
account in the present work. However, within the opti- 
mally doped to overdoped regime, the superfluid density 
is large and the critical temperature corrected by these 
latter fluctuations is expected to be close to the mean- 
field BCS one, even when the transition is Kosterlitz- 
Thouless rather than BCS-like jl5). For dopings accessi- 
ble experimentally in the cuprates (x < 0.5) the s-wave 
superconducting instability is suppressed. Including the 
effect of antiferromagnetic spin fluctuations (AFSF), we 
find that three different regimes remain as a function 
of doping in the realistic range: (i) long-range antiferro- 
magnetism near half filling, (ii) d-wave superconductivity 
in a pure Fermi liquid (usual BCS) scenario within the 
overdoped regime; (iii) d-wave superconductivity in the 
presence of AFSF within the underdoped to optimally 
doped regime. The maximum T c within our BCS treat- 
ment depends on the strength of the interaction and the 
position of the VHS. The effect of the spin- fluctuations, 
however, modifies the BCS picture within this regime. 
The effective pairing interaction of (Q) has the form of 
an exchange coupling, being ineffective in the Neel state 
p6f , but providing an explicit coupling in the regime of 
doping for which the AF correlations are short ranged. 
We considered particular values of the model parame- 
ters for which optimal doping, T c and shape of the Fermi 
surface are in agreement with experiments and further 
examined the effect of the non-trivial temperature be- 
havior induced by the existence of hot spots in the Fermi 
surface in the presence of AFSF. We found that T c de- 
creases as the correlation length increases in correspon- 
dence with experiments. For optimally doped materials, 
we expect that the magnitude of T c is dominated by the 
effect of spin fluctuations rather than by phase fluctu- 
ations of the superconducting order parameter. Details 
of the coupling with AFSF, renormalizations of the bare 
interactions, and eventual consequences upon the pseu- 
dogap behavior requires a treatment of the full T-matrix 
of the effective pairing interaction on equal footing with 
the spin susceptibility and is left for future studies. 
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